Nonmonotonic charge occupation in double dots 
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We study the occupation of two electrostatically-coupled single-level quantum dots with spinless 
electrons as a function of gate voltage. While the total occupation of the double-dot system varies 
monotonically with gate voltage, we predict that the competition between tunneling and Coulomb 
interaction can give rise to a nonmonotonic filling of the individual quantum dots. This non- 
monotonicity is a signature of the correlated nature of the many-body wavefunction in the reduced 
Hilbert space of the dots. We identify two mechanisms for this nonmonotonic behavior, which 
are associated with changes in the spectral weights and the positions, respectively, of the excitation 
spectra of the individual quantum dots. An experimental setup to test these predictions is proposed. 
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Introduction. - Quantum dots coupled to an elec- 
tron reservoir can be filled with electrons one by one 
in a controlled way. With increasing gate voltage, the 
number of electrons on the dot increases in a step-like 
manner. As long as Coulomb interaction among the dot 
electrons is negligible, the dot levels are filled indepen- 
dently from each other. However, small semiconductor 
quantum dots are usually subject to strong Coulomb in- 
teraction: charging energy gives rise to extended plateaus 
of (almost) fixed dot charge as a function of gate volt- 
age. The aim of the present paper is to demonstrate that 
more dramatic signatures of Coulomb correlations can 
be found in studying the filling of the individual levels 
in the quantum dot. While the total occupation of the 
dot remains a monotonic function of the gate voltage, we 
predict that, under circumstances specified below, the 
individual levels are filled in a nonmonotonic way. The 
physics proposed here enables one to observe non-trivial 
many-body correlations in a particularly simple (effec- 
tive) 4-dimensional Hilbert space. 

To illustrate the signature and mechanisms of the non- 
monotonic filling, we consider the simplest possible model 
system which allows for a separate access to the individ- 
ual level occupation: two electrostatically-coupled single- 
level quantum dots with spinless electrons, see Fig. ^ 
The double dot is equivalent to a single quantum dot ac- 
commodating two levels. The advantage of the double- 
dot setup lies in the possibility to read out separately 
each dot's occupation by electrostatically-coupled quan- 
tum point contacts. This kind of charge sensing has been 
experimentally demonstrated for single 1] and double 
dots recently. The quantum dots are tunnel coupled 
to one common or to two separate electron reservoirs. 

We model the double-dot system with the standard 
tunnel Hamiltonian H = H dot i + H dot 2 + H ch + flieadi + 
-Hiead2 + Ht,i + Ht,2- Quantum dot i = 1,2, de- 
scribed by Hdot,i — £ic\ci, accommodates a level with 
energy ej, measured relative to the Fermi energy in the 
leads. The level energies can be tuned by the gate 
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FIG. 1: Two electrostatically-coupled single-level quantum 
dots are attached to electron reservoirs and gated by gate 
electrodes. The individual occupation is determined by adja- 
cent quantum point contacts. 



electrodes. Without loss of generality, we always as- 
sume £2 > ei. The charging energy is accounted for by 
-Hch = Un\ri2- Each dot is coupled to an electron reser- 
voir 0|, -HjeacU = J2k ewo^afcj. Tunneling is modeled by 
i?T,i = J2k {^i c l a ki + h.c.J , where we assume the tunnel 
matrix elements to be energy independent. The tunnel 
coupling introduces a line width Ti = 2n\ti\ pi, where pi 
is the density of states of lead i at the Fermi energy. 

There are four possible states for the double-dot sys- 
tem: the double dot being empty (x = 0), singly (\ — 
1,2), or doubly (\ = d) occupied. The corresponding 
energies are E x [4( . To address the question of how the 
occupations (ni) and (n 2 ) of the dots vary as both levels 
ei and 62 are simultaneously pulled down, we keep the 
bare level separation A = £2 — e% fixed. The fact that tun- 
neling does not commute with the kinetic and charging 
terms is essential for the physics outlined below. 

Regimes. - We begin with identifying the regime at 
which a nonmonotonic filling of the individual dots are 
expected. Consider first [7 = 0. The spectral density 
of dot i is a Lorentzian centered around e» with a width 
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IV As ti crosses the Fermi energy of the lead, the oc- 
cupation (rii) of level i changes from to 1. The width 
of this transition is governed by max-jTj, k^T}. It is ob- 
vious that in this case the occupation of each level is a 
monotonic function of the gate voltage. We, thus, turn 
our attention to the limit U 3> max{A, V, k-gT}. 

For temperatures larger than the level broadening, 
fceT ^> T, the influence of the tunnel couplings on the oc- 
cupations is negligible and the double-dot system is filled 
according to the Boltzmann factors exp(— £7^ /feeT) of 



the corresponding bare energies E x ^ with ej = Ef'—E^ 
and U = Ef ] -Ef ] -Ef ] +Ef\ Now, the occupations 
of the two dots are correlated due to charging energy, 
(nin 2 ) — (ni)(n 2 ) < 0, but the individual dot occupa- 
tions (rii) are still found to depend monotonically on the 
gate voltage. Therefore, we require T > k^T. 

If the level separation is small as compared to either 
level width or temperature, max{fcBT, T} ^> A, then the 
double dot will be filled in a symmetric way, (n%) w (71-2}, 
which again yields a monotonic filling. This restricts the 
regime of interest to A > T > k-gT. 
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FIG. 2: For k B T > F or F » A, the occupation of the 
individual dots is a monotonic function of the gate voltage. 
Nonmonotonicities due to mechanism 1 are strongest for A ~ 
r 2> fceT, and mechanism 2 becomes effective for A ~ V ~ 
k B T and T 1 ^T 2 . 

Mechanisms for nonmonotonicity. - We identify 
two different mechanisms by which a nonmonotonic fill- 
ing of the individual levels can occur. They are both 
related to the Coulomb interaction, with two different 
consequences concerning the spectral density. First, as 
mentioned above, there is spectral weight of level 2 not 
only around £2 (describing filling or depleting of dot 2 
while dot 1 is empty) but also around £2 + U (describing 
transitions when dot 1 is filled) . The total spectral weight 
is normalized to 1, with the relative spectral weight be- 
tween the two peaks being determined by occupation of 
dot 1. Both peaks of level 2 are broadened by T2 |H, 
which yields a partial filling of the dot once the peak is 
close to the Fermi energy. As the level in dot 1 passes 



through the Fermi level and dot 1 is filled, spectral weight 
for the level in dot 2 is transfered from the peak around 
£2 to the one around £2 + U, which is further away from 
the Fermi level, yielding a reduction of the partial oc- 
cupation of dot 2, giving rise to a nonmonotonic filling. 
This mechanism will be visible for A ^> kBT (since oth- 
erwise the filling is dominated by thermal fluctuations) 
and is most pronounced for A ~ T ^> kBT. 

The second mechanism is based on a renormalization of 
the peak positions in the spectral density due to Coulomb 
interaction 0. As discussed in more detail below, tun- 
neling in and out of dot 2 yields a renormalization of 
dot level 1, t\ — > e±. The renormalization is strongest 
when level 2 is close to the Fermi energy. This renormal- 
ization can give rise to a nonmonotonic filling of dot 2 
under the condition A ~ fceT and (rii) + (71,2) ~ 1. In 
this case, both dots will be partially filled once £1 and 
£2 are below the Fermi level. The relative occupation 
is approximately given by exp(— A/fceT). A renormal- 
ization of the level splitting A then changes the relative 
occupation. If level 2 is close to the Fermi level, then 
level 1 is strongly renormalized towards higher energies. 
As a consequence, the occupation of dot 2 is increased 
on the cost of level 1, i.e., there is a reshuffling of the 
relative occupation. When the gate voltage is increased 
further, the renormalization becomes weaker, and level 
1 gains back some of the occupation it lost to level 2, 
which leads to a nonmonotonicity of (n 2 ). It turns out 
that for symmetric coupling, Ti = T2 = T, the renormal- 
ization of the level separation is negligible in the regime 
A ~ r ~ kBT since both levels are renormalized approx- 
imately by the same amount. To get a nonmonotonic dot 
filling, an asymmetric coupling Ti 7^ T2 is required. 

Perturbation expansion. - To substantiate the 
qualitative ideas outlined above, we perform a pertur- 
bation expansion of the dot occupations in the tunnel- 
coupling strength T. For this, we employ a diagrammatic 
imaginary-time technique that was introduced to study 
charge fluctuations in a metallic single-electron box [3, 
adjusted to a two single-level-dot system. The central 
quantity is the partition function Z = ^2 x cxp(—[]E x ). 

The corresponding effective energies E x = E x ^ + E x ^ + 
C(r 2 ) can be expanded in orders of T. Evaluation 



of E x 1] yields E. 
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-Mei) + <7 2 (e 2 + U)}, 4 1} = -[oi(ci + U) + <7 2 (e 2 )], 
and Ef ] = -[oi(ei + U) + a 2 {t 2 + U)], with a t {uj) = 
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where the bandwidth 



D appears as a high-energy cutoff at this intermediate 
step but drops out for all physical observables since they 
only depend on differences of E x 's. As a result, the dot 
levels are renormalized according to £; = £; 
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for the same expression holds with 1 replaced by 2, 
and U remains unrenormalized in first order in T. Here, 
^>(x) is the digamma function. From Eq. JJ) we see that 
the renormalization of the level in dot 1 is proportional 
to the tunnel coupling of dot 2 (and vice versa). It is 
maximal when e 2 or e 2 + U is in resonance with the Fermi 
level, and vanishes in the absence of Coulomb charging. 

Thermodynamic quantities such as the occupation of 
the dots are obtained by adequate logarithmic derivatives 
of the partition function. In particular, 
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The first-order corrections of these quantities are directly 
obtained from 
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l-(«2) (0) Ja 1 (e 1 ) + (n 2 )(°) ( 7 1 (e 1 + [/) 
1 - <m) (0) ] <7 2 (e 2 ) + <«i) (0 V 2 (e 2 + U)] (4) 



plus a constant that is independent of ei, e 2 , and U. 

The set of Eqs. ©-© provides the starting point for 

the subsequent quantitative analysis. For example, the 
correction to occupation of dot 2 is 
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where we have used d(n 2 )^ /dei = d{ni)^ / de 2 . 

Results. - In Fig.|3|we show the dot occupations as 
a function of the mean level position e = (ei + e 2 )/2 in 
the regime A ~ T » k B T. We find nonmonotonicities 
which we relate to mechanism 1 discussed above. For 
an approximate analytical understanding of the result, 
we concentrate on the nonmonotonicity of (n 2 ) on the 
right step. Since e 2 > ei and A » k B T, we can set 
(n 2 )(°) « for the whole region over which level 1 is 
filled up, and (ni)^ is independent of e 2 . Nevertheless, 
level 2 is partially filled due quantum fluctuations, 
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As level 1 passes through the Fermi level, filling up quan- 
tum dot 1, the occupation of dot 2 drops from r 2 /(27re 2 ) 
down to r 2 /[27r(e 2 + U)]. Within our perturbative anal- 
ysis, the width of this drop is provided by temperature, 
as a broadening of level 1 due to quantum fluctuations 
would enter in higher orders only. When going to higher 



order in T, the width is expected to be max{A;BT, Ti}. 
The amplitude of the drop is approximately 
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FIG. 3: Nonmonotonic dot filling due to mechanism 1. The 
parameters are A = 0.1(7, I\ = T 2 = 0.1(7, and k B T = 0.01(7 
[3- Inset: Result for the symmetrized Hartree approximation. 

We now turn to the regime A ~ T ~ k-g,T to dis- 
cuss mechanism 2. From inspection of Eq. Q we real- 
ize that in the given regime the renormalization of both 
levels is roughly the same. In order to generate a sub- 
stantial renormalization of the level separation that can 
give rise to a nonmonotonicity, we choose different cou- 
pling strengths Ti ^ T 2 . In the following we discuss 
the case T 2 > T\. The result is shown in Fig. 01 The 




FIG. 4: Nonmonotonic dot filling due to mechanism 2 for 
A = 0.117, Yi = 0.01(7, T 2 = 0.1(7, and k B T = 0.0517. 

second line of Eq. (JSJ) is associated with the renormal- 
ization of level 1. The renormalization (towards higher 
energies) is strongest when e 2 passes through the Fermi 
energy. The peak of the nonmonotonicity is, therefore, 
located to the left as compared to the one related with 
mechanism 1. The width is given by T 2 . For an esti- 
mate of its height S(n 2 ), we neglect both the possibil- 
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ity of the double dot being empty or doubly occupied, 
i.e., {712)^ ~ exp(— /3A)/[exp(— /3A) + 1], and we assume 
U > k B T. We find 

^ " ( e -gA + 1)2 2 ^ ( 8 ) 

If we reverse the asymmetry of the coupling strengths, 
Ti > T2, we find, by using analogous arguments as above, 
a nonmonotonicity of (ni) on the left step, while (712) 
remains monotonic. 

In Fig. [3] we show the correlator (niria) — {n\){n2), 
expanded up to first order in T, for the two set of pa- 
rameters chosen in Fig. [21 and Fig. 0] We find that the 
nonmonotonicity of the dot occupation is accompanied 
with an enhancement of correlation. 
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FIG. 5: Correlator (711.712) — (ni)(ri2). The parameters are as 
in Fig. [3] (solid line) and Fig. [I] (dashed line). 

Throughout our analysis, we assumed flat bands, i.e., 
particle-hole-symmetric densities of states, in the leads. 
As a consequence, (ni) at e for given r^Ta is identical 
to 1 — (712) at U — e with the values of T\ and T 2 being 
interchanged. Thus, our discussion, which was focused 
on the right step of (712) holds for the left step of (n\) as 
well (after interchanging T\ T2). 

In parallel to this perturbation theory, we have per- 
formed an equation-of-motion (EOM) analysis up to the 
second hierarchy solving a set of 1 1 coupled equations for 
the various Green's functions. This approach, however, is 
not systematically controllable. Truncating the EOMs on 
the first hierarchy results in the symmetrized Hartree ap- 
proximation with the Green's function G t ii(oj) given 
by [Ggiw)]- 1 = [(1 - (n a » /(w - ei) + (na>/(w - e x - 
U)] 1 +iT\/2 and similar for G^(iJ). The resulting (m) 
and (712) tend to be biased towards each other, but the 
nonmonotonicity is still observed, see the inset of Fig. [3] 

Summary. - In summary, we predict that the com- 
petition between tunneling and charging energy can give 
rise to nonmonotonic filling of individual quantum-dot 
levels. We identify two different mechanisms leading to 
nonmonotonicities and determine the regimes at which 
they occur. Based on a perturbation expansion in the 



tunnel couplings, we derive analytic expressions for the 
location and the strength of the signal. 
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